AD-A268  111 


NAWCWPNS  TP  8098 


An  Introduction  to  Wideband  Ambiguity 
Functions  and  Time-Scale  Representations 
for  Radar/Sonar  Applications 


by 

Brett  Borden 
Physics  Division 
Research  Department 


JANUARY  1993 


NAVAL  AIR  WARFARE  CENTER  WEAPONS  DIVISION 
CHINA  LAKE,  CA  93555-6001 


9 


Approvsd  for  public  release;  distriHition  is  unimiied. 


93-20081 


Naval  Air  Warfare  Center  Weapons  Division 


FOREWORD 


The  research  described  in  this  report  was  poformed  during  fiscal  year  1993  as  part  of 
an  effort  to  in^rove  upon  radar  classification  and  identification  of  noncooperative  targets. 
This  problem  continues  to  be  a  primary  goal  of  radar  research  programs,  and  considerable 
effort  has  been  expended  witUn  the  l^t  few  decades  in  attempts  to  solve  it  Invariably, 
the  accuracy  of  these  methods  depends  upon  available  radar  resolution  and  so,  indirectly, 
it  depends  on  radar  bandwiddi. 

Many  classical  results  in  radar  systems  theory  presume  that  the  narrowband 
approximation  is  valid.  It  is  hop^  that  dtis  r^ort  wiU  help  clarify  what  is  gained  and 
what  is  lost  when  this  approximation  is  not  valid.  The  results,  while  not  surprising,  also 
offer  an  interesting  interpretation  in  terms  of  "wavelet”  (or  time-scale)  representations. 
This  effort  was  supported  by  the  Office  of  Naval  Research. 

This  report  has  been  reviewed  for  technical  accuracy  by  Gary  A.  Hewer. 


Approved  by 

R. L.  DERR,  Head 
Research  Department 
15  January  1993 

Released  for  publication  by 

S. HAALAND 

Deputy  Commander  for  Research  &  Development 


Under  authority  of 
W.  E.  NEWMAN 
RAdm.,  U.S.  Navy 
Commander 


NAWCWPNS  Technical  PubUcation  8098 


Published  by .  Technical  Information  Department 

Collation .  Cover,  17  leaves 

First  printing . 65  copies 


REPORT  DOCUMENTATION  PAGE 


Fdm  Approved 
OMB  No.  0704^188 


Pubic  wpmtiwgbufdtn  far  llii»oollialione<in(oimNian««NimNiJ  to  awwAga  1  hour  per  l■ipon»•.(■olUdinBlh•lim*tan•»i■*ingil■(nl^lion».^••fc^linB•«■liBgda••ou»o»,#N^l•^iB9«>d 
mcinloininqlhodaomodod.oiidootnilolingiidiouiwiwiglhoooloehonolinloniMlion.  Bond  oommonti  ragoiding  thio  buibon  oaimalo  or  any  other  apace  ci  Mo  oolcdion  d  Monmbon, 
inokidjns  luggoNiono  for  raduoing  Mi  budon,  to  Waohington  Hoodqiiitoii  Sorvino,  DNotonto  lor  kdonmion  Opmliont  and  Ropodi.  1216  Mloraon  OoMt  HighuMy,  Suio  120*.  Aiington, 

VA  22202-4302,  told  to  thoOflioaai  Maragonwil  told  Budgal.Papo>woifcRodiiolionPraiHt(070«-01M),Wiahiti^,  DC  20603. 

1.  AOENCV  USE  OiLV  (lm¥»  blank) 

2.  REPORT  DATE 

January  1993 

3.  REPORT  TYPE  AND  OATES  COVERED 

Interim  Oct— Dec  1992 

4.  TITLE  AND  SUITITLE 

AN  INTRODUCTION  TO  WIDEBAND  AMBIGUITY  FUNCTIONS  AND 
HME-SCALE  REPRESENTATIONS  FOR  RADAR/SONAR  APPLICATIONS 

&  FUNONM  NUMBERS 

WU 113132 

PE  61153N 

TA  BR052-Q2-01 

*.  AUTHOIKS) 

Brett  Borden 

7.  PERFORMNO  OROAMZATION  NAME(S)  AND  A0DRE8S(ES) 

Naval  Air  Warfare  Center  Weapons  Division 

China  Lake,  CA  93555-6001 

IL  PERFORMINQ  ORGANIZATION 

REPORT  HUMER 

NAWCWPNS  TP  8098 

b.  SPONSORINO/HdUTMlUfl  AOENCY  NAME(S)  AND  ADBRESS<E8) 

Office  of  Naval  Research 

Arlington,  VA  22217-5000 

NL  SPONSOflINQ/MONITORINQ 

AGENCY  R^ORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

1  ItA.  nSTTHBUTION/AVAILAllLnY  STATEMENT 

Its.  DISTRIBUTION  CODE  I 

A  Statement;  public  release;  distribution  unlimited 

1  13.  ABSTRACT  (Muamum  200  woidi) 

(U)  Recent  advancements  in  the  ability  of  radar/sonar  ^sterns  to  laundi  and  analyze  signals  of  very  short 
time  duration  have  caused  many  to  re-examine  the  so-called  wideband  ambiguity  function.  Much  of  this 
work  has  been  at  a  mathematical  level  that  has  placed  it  out  of  reach  to  many  radar  and  sonar  engineers. 

Some  of  the  work  has  been  incomplete— occasionally  even  incorrect  Our  present  pedagogical  discussion  uses 
the  standard  time-frequency  motivation  for  time-scale  distributions.  We  compare  and  contrast  the  usual 
narrowband  techniques  wiffi  the  wideband  approach  and  explain  how  high-resolution  downrange  imaging 
requires  an  understanding  of  the  latter. 

14.  SUBJECT  TERMS 

Ambiguity  function 

16.  NUMBER  OF  PAGES 

32 

wiQeD&nc 

Wavelet 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

16.  SECURITY  CLASSinCATION 

OF  THM  PAOE 

UNCLASSIFIED 

It.  SECURITY  CLASSIRCAT10N 

OF  ABSnUCT 

UNCLASSIFIED 

ta  UMTATION  OF  ABSTRACT 

U. 

Standard  Fonn  298  (Rav.  2-89) 

PiMOriM  by  ANSI  SW.  23S-U 

2M-ioe 


NSN  7540-01-280-5500 


NAWCWPNS1P8098 


CONTENTS 


1.  Introduction .  3 

2.  Problem  Statement .  3 

3.  Derivation  of  P(t,<D) .  4 

4.  Kernel  Function .  7 

Example  4.1 .  8 

Exanq;>le  4.2 .  9 

Example  4.3 .  9 

Exan^le  4.4 .  9 

5.  (Radar)  Cross-AmhiguiQr  Function .  12 

6.  Choice  of  Signals .  16 

7.  Why  Go  Wideband? .  24 

8.  Conclusion  .  27 

Appendix 

Motivation  for  the  Measure  dT  dDID .  29 

References .  31 


Ij  :  tc  T/AI-XTY  irrSfECTEI;  3 


Aataaalen  For 


WtlS  ORAftl 
DIXC  TAB 
VDaan0UD0»a 
JttStlf l«8ll0fU 


□ 

□ 


By 

Availabllltf  flsyXM 

iTBll 

^0 

Bad/or 

Lai 

NAWCWPNSTP8098 


1.  INTRODUCTION 


Within  the  last  few  years  there  has  been  a  renewal  of  interest  in  wideband  radar 
systems  and  techniques  (References  1  tluough  3).  This  interest  is  largely  because  of  the 
prospects  for  high-resolution  downrange  target  imaging  that,  because  of  technological 
advancentents  in  signal  generation  and  processing,  are  o^y  now  becoming  feasible.  It  is 
also  due  to  inherent  differences  between  the  natural  representations  of  narrowband  and 
wideband  radUr  signals  and  a  new  appreciation  of  how  these  differences  might  be 
exploited. 

The  purpose  of  the  present  discussion  is  to  examine  some  of  the  results  from  this 
field  that  are  relevant  to  the  applied  problem  of  increased  resolution  in  radar  target 
imaging.  Much  of  the  recent  woik  in  wideband  ambiguity  frmctions  has  had  a  decidedly 
maAematical  flavor  that,  regrettably,  has  put  it  out  of  Ae  reach  of  many  of  the  radar 
engineers  who  might  be  expected  to  tenefit  most  from  it  We  shall  avoid  this  iq>proach 
whenever  possible  and  choose  explanations  that  are  more  "classical"  in  appearance.  This 
means  that  we  will  invariably  be  less  accurare  in  our  definitions,  and  the  reader  is  advised 
to  refer  back  to  the  original  references  for  coirpleteness. 

We  begin  by  briefly  reviewing  the  standard  time-frequency  methods  in  signal 
processing.  We  shall  m^e  no  attempt  to  reference  all  of  the  important  achievements 
made  by  Ae  many  researchers  in  this  field;  instead,  we  refer  the  interested  reado*  to  the 
excellent  review  of  this  material  presented  in  Reference  4.  In  Section  S  we  introduce  the 
wideband  and  narrowband  ambiguity  functions.  These  functions  are  keys  to 
understanding  radar  imaging,  and  we  will  show  how  they  are  related  and  how  they  differ. 
While  the  narrowband  ambiguity  function  is  understandable  in  terms  of  time-fr^uency 
methods,  we  will  see  that  the  wideband  ambiguity  function  requires  time-scale 
techniques.  In  Sections  6  and  7  we  will  explore  some  of  the  differences  between  time- 
frequency  and  time-scale  representations  and  argue  for  choosing  the  latter  in  high- 
resolution  imaging  applications. 


2.  PROBLEM  STATEMENT 


For  a  signal  s{t),  denote  by  \s(tj^  the  "instantaneous  power"  at  time  t  and  by  jiC  (oj^ 

the  "power  spectrum"  at  ca.  We  seek  to  determine  a  function  called  the  time- 

frequency  distribution,  which  satisfies  the  marginal  relations 
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«•  «• 

J  P(t,o})do) = Isrrf .  J  P(t.to)dt = l57o)j|* . 


(2.1) 


and  such  that  the  total  eneigy  of  the  signal  is  given 


£=  J  jP(t,ca)do)dt . 


(2.2) 


For  stochastic  signals,  these  results  reduce  to  their  probabilistic  form.  For 
deterministic  signals,  this  distribution  is  really  an  energy  density;  the  "distribution" 
nomenclature,  though  partly  a  holdover  from  quantum  mechanics,  is  used  because  the 

global  average  (the  expectation)  of  a  function  g(t.n>)  is  defined  by 


(«>=  J  jg(t.(o)P(t,(o)dtd(0  . 


(2.3) 


The  local  or  conditional  mean  at  a  particular  time  is  defined  by 

«• 

jg{t,Q))P{t,a})d0 


(and  similarly  for  the  mean  at  a  particular  frequency). 

Equations  2.3  and  2.4  are  the  main  motivations  for  seeking  P{t,0). 


(2.4) 


3.  DERIVATION  OF/>(r,Q)) 


f 

The  characteristic  function  is  the  expectation  ofe'^****^; 


M{0.t)s  J  P{t.0)dtd0  . 


(3.1) 
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Characteristic  functions  are  convenient  because  joint  moments  can  be  calculated  by 
differentiation: 


<<■<»■)  =  7= 


.t«0 


(3^) 


and  so,  expanding  the  exponential  in  Equation  3.1  in  a  Maclauiin  series  reveals  that  they 
may  be  expressed  as 


hbOmvO 


(3.3) 


and 


For  functions  g,(r)  and  £2(0)),  we  have 

ta(0)=  J  J  s*(e))g,^i-^^s(o))do) 


(3.4) 

(3.5) 


Consequently,  we  can  associate  time  and  frequency  with  die  opoators  T  and  W ,  so  that 


r  ->  r  IV  -I—  in  the  time  domain 
dt 


T->i-^  W^O) 
dCD 


(3.6) 


in  the  frequency  donuiin. 


A  function  g(t,0))  of  both  time  and  frequency  may  be  formally  treated  in  the  same 


way: 


(g(r.fl)))=  ]s*(0G(My)s(/)dr=  Ji*(o>)G(r.ai)f(a>)dn>  , 


(3.7) 


where  G{T,W)  is  the  operator  "associated"  with  g{t,(o).  Equation  3,3  allows  us  to  write 
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M{0.  T)  *  ->  J s\t)e^  e‘^  $(t)  dt 

or  (3.8) 

However,  there  is  an  inherent  ambiguity  in  Equation  3.8  because  the  association 
g(t.ct))  G{T,W)  is  not  unique.  (For  example.  Equation  3.3  also  allows  us  to  choose 
^•(•f+w)  ^.(ar+tw)  ^  unambiguous  procedure  for  the  association  of  G{T.W)  with 

g{t.  Of)  sets  (References  4  and  5) 


G(T.W)=J  ,  (3.9) 

where 

y(®.  t)  =  ^  J  J g(t.  dtd0  (3.10) 


and  ^(0,  t)  is  a  function,  called  the  kernel,  that  must  satisfy  certain  initial  conditions  if 
the  correct  marginals  are  to  be  obtained  (see  below). 

In  the  time  domain,  the  operator  is  the  translation  operator  [i.e., 

e‘^  s(t)  =  s(t)  -  s(f  +  t)]  and  substitution  results  in  the  general  form  for  the 
characteristic  function  as 


M(0,t)  =  ^(0,T)js*(u)e’*‘s(u+T)du 


(3.11) 


Inverting  ^nation  3.1  with  M(0,  t)  given  by  Equation  3.1 1  yields  the  associated  genoal 
form  for  distributions: 


^(0,  t)s*(u)s(u + x)dudzd0. 


(3.12) 
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There  are  other  forms  that  are  used  finequently  to  define  this  so-called  Cohen's  class 
of  tiine-fiequency  distributions.  For  exanq>le,  we  can  write 


F(r,n))  =  ^  J  ]w{u.o))(l>{e.T)e-**^-'‘^dude  ,  (3.13) 

where 

W(u,a))  s  ?  s*(u)s(u  +  dr  (3.14) 


is  the  Wigner-Ville  distribution.  In  this  form,  the  factor  0(0,  may  be 

interpreted  as  a  "smoothing  taper." 

Another  form  writes 

F(/.fi))  =  ^  j  J;f(0.T)0(0,T)e“'^**”*>d0dT .  (3.15) 


where 


X{^>  t)  s  ?  5*(«)s(«  +  r)e'^  du 


(3.16) 


is  the  (narrowband)  ambiguity  function,  which  will  be  discussed  in  detail  below. 


4.  KERNEL  FUNCTION 


The  propCTties  of  the  distribution  function  are  a  direct  consequence  of  the  choice  of 
kernel.  This  kernel,  in  turn,  cannot  be  chosen  at  will  but  must  obey  certain  integration 
and  transformation  criteria.  In  particular,  the  kernel  determines  whether  the  density  is 
correctly  related  to  the  instantaneous  energy  and  spectrum  (References  4  and  5). 

Integration  of  P{t.a>)  with  respect  to  0  yields 
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J  P{t,o)tUa  -  t)s‘iu)s(u  +  i)d0dudT 

If  this  is  to  equal  .  then 

±  ^{e,G)de  =  8{t-u)  =»  0(0,0)  =  1 . 

2X^ 


(4.1) 


(4.2) 


Similaily,  it  is  easy  to  show  that  0(0,  t)  « 1. 


We  assume  throughout  that  5(t)  is  normalized  to  unit  energy,  and  so  it  also  follows 
that 

J  P{t,(o)daidt  *  1  =  total  enagy  =*  0(0,0) » 1 .  (4.3) 


If  we  require  P{t,(o)  to  be  teal  (as  is  appropriate  for  a  description  of  energy 
density),  then  0(0, t)  =  0*(-0,-t). 


EXAMPLE  4.1 

If  we  choose  0(0,  r)  =  1,  then  we  obtain  the  Wigncr-Ville  distribution  (References  4, 
6,  and  7); 


P^{t,(o)  =  ^  \e-^^s'‘(t)s(t+  T)dt .  (4.4) 

The  Wigner-Ville  distribution  has  been  very  well  studied,  and  it  is  known  that  it  may 
introduce  reconstruction  artifacts  for  multicon^nent  signals  (References  4, 8,  and  9).  In 

particular,  this  distribution  suffers  in  that  it  may  not  be  zero  even  when  the  signal  is.  ♦ 
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EXAMPLE  4/* 

The  Choi-Williams  distribution  results  from  setting  0(0,  t)  =  «"**** "  so  that 
(References  4  and  9) 

J  J  +  T)dudT  .  (4.5) 

The  Choi-WUliams  distribution  can  be  thought  of  as  a  "smoothed"  Mgner- Ville  dis¬ 
tribution  with  the  parameter  tr  chosen  to  control  spurious  cross-terms.  ♦ 


EXAMPLE  4  J 

The  choice  0(0.  t)  =  e***  results  in  the  Kiriewood-Rihaezek  distribution  (References 
4, 10,  and  1 1); 

.  (4.6) 


This  is  sometimes  called  the  complex  energy  spectrum.  ♦ 

The  distributions  of  Examples  4.1  through  4.3  result  in  the  correct  marginal 
distributions.  Moreover,  the  Wigner-Ville  and  Qtoi- Williams  distributions  are  also  real 
valued.  However,  these  (Ustributions  generally  are  not  positive  and  it  may  be  difficult  to 
interpret  expectations  based  upon  them  (Reference  4).  A  distribution  that  is  always 
positive  is  the  spectrogram. 


EXAMPLE  4.4 
The  spectrogram. 


^sp«.(^Q>) 


MT)h(T 


2 

-t)dx  , 


is  obtained  fit>m 

•• 

0(0,  t)  =  J e-**  h{u  +  x)h*{u)du  . 


(4.7) 


(4.8) 
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The  spectrogram  is  generally  an  easy-to-inteipret  distribution.  However,  its  flaw  is 
that  it  does  not  yield  the  correct  marginals  (satisfying  ^(0.0)  =  0(0,  r)  =  1  requires  that 

|A(f)f  =  5(r)  and  |A(a))j  =  5(fl)),  and  this  is  impossible).  In  practice,  the  "window"  h{t) 

is  tailored  to  a  particular  application  with  looking  "more  like"  a  delta  function  if 

good  time  resolution  is  desired  and  |A((ii})|  looking  more  like  a  delta  function  if  good 

frequency  resolution  is  needed.  ♦ 

A  positive  distribution  that  does  yield  the  correct  marginals  is 

/»(/,©)  =  |s(f)f|s(Q))f.  (4.9) 


This  distribution  is  determined  only  by  its  marginals  and  so  lacks  time-riequency 
correlation  information  (resulting  in  possible  artifacts  in  signal  reconstruction).  In 
addition,  this  distribution  is  not  bilinear.  In  fact,  it  is  impossible  to  have  positive  joint 
distributions  with  correct  marginals,  which  are  bilinear  in  s(t).  However,  it  is  possible  to 
satisfy  positivity  with  coirect  marginals  when  0(0,  t)  is  a  functional  of  j(f).  In  this  case. 
Equation  4.9  can  be  generalized  (References  12  and  13),  although  a  systematic  procedure 
for  incorporating  the  correlation  information  has  not  been  developed  Reference  4). 

For  any  d^ired  expectation,  the  spectrogram  will  yield  only  an  average  over  the 
extent  of  the  window  and  can  be  relied  upon  oiidy  when  the  signal  is  short-time  stationary 
(Reference  8).  While  the  s|^trogram  does  not  satisfy  the  marginals,  its  positivity  is  a 
very  useful  feature.  In  addition,  it  is  known  that  the  derivative  of  the  phase  does  not 
dways  correspond  to  the  frequency  in  the  Fourier  spectrum  and  so  there  may  be  an 
inherent  inconsistency  with  our  original  requirement  that  the  distribution  satisfy  the 
margin^  constraiiits  References  4  and  14).  Consequently,  in  the  following,  we  shall 
prefer  interpretability  over  "correctness" — especially  since  it  is  recogniz^  that  the 
current  theory  is  not  fully  consistent  and  requires  considerable  further  research 
(Reference  4). 

The  behavior  of  the  kernel  and  associated  distribution  functions  under  various 
transformations  of  the  signal  are  listed  in  Table  4.1  (c.f..  Reference  4). 
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TABLE  4.1.  Properties  of  Distributions  for  Transformations  of  the  Signal. 


s(t) 

lime  shift 

P(t+t^,a)) 

Any 

Frequency  shift 

s{a) + u), )  or  f(r)e'**' 

P(t.Q)  +  0)„) 

Any 

Hme  scaling 

^|\^s{oa) 

P(ea,o)ta) 

<ft{e  1  a.at)  = 

Freq.  scaling 

P(i/p,pa) 

Hme  inversion 

5(-r) 

P{-t,-(o) 

0(-0,-t)  =  ^(B,t) 

Complex  conj. 

P{t.-(0) 

Finally,  and  for  completeness,  we  need  to  discuss  the  so-called  uncertainty 
principle.  This  relationship  between  the  standard  deviation  of  a  function  and  the  standard 
deviation  of  its  Fourier  transform  is  (Reference  4) 


V— 


2 


(4.10) 


for  any  signal  s(t). 

The  uncertainty  principle  is  a  relationship  concerning  the  marginals  and  expresses 
the  well-known  result  that  a  function  and  its  Fourier  transform  cannot  be  made  arbitrarily 
narrow. 
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5.  (RADAR)  CROSS-AMBIGUITY  FUNCTION 


For  a  point  target  at  distance  r  from  a  stationary  radar  and  moving  with  downrange 
speed  V,  Iv  I «  c,  the  echo  returned  from  a  transmitt^  pulse  w(t)  is  given  by  (R^erences 
1  through  3,  IS  through  17) 


(5.1) 


(The  scale  factor  is  required  for  conservation  of  energy.)  This  is  an  exact  result  and  can 
be  written  as 


e(r)  =  VDH<D(r+r)) , 


(5.2) 


where  D  =  (c-v)l(c+v)  is  the  "Doppler  stretch  factor"  and  T  =  -2rl{c-v)  is  the 
signal  delay  at  r  =  0. 

The  cross-correlation  of  a  test  pulse/,  =Vflw^(a(«+r))  with  a  point  target  echo 
e^{t)  from  w,(0  is 


(«,./«)  =  JVd  w,(D(«+r))Va  w;(fl(«+f))dM 


(5.3) 


If  we  defrne  the  wideband  cross-ambiguity  function  by  (References  1  through  3, 17) 
T,a)  s  Va  J w^(a(u  +  t))w*  (ii)dM ,  (5.4) 

then  we  can  write 

(e../.Xt.a)  =  A_(<i(r  -  »),f )  .  (5.5) 
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The  narrowband  cross-ambiguity  function  ^  obtained  from  Equation 

5.4  by  writing  both  w^  and  w,  in  the  foim  w{t)  =  where  0)„  may  be  inteipreted 

as  a  carrier  frequency  ^ut  need  not  be).  Substitution  yields 


«■ 

Amir, a)  =  Vae'*'"  j s,(a(u+  r))sl(u)e‘^du  , 


(5.6) 


where  0  s  (a  - 1)©,.  Expanding  f,(a(M+  t))  in  powers  of  e  =  2v  /  (c + v)  and  dropping 

all  but  the  lowest  order  term  (i.e.,  restricting  s{t)  to  be  slowly  varying  with  t  in 
comparison  to  the  exponential  factor)  results  in  (Reference  17) 


[Note  that  the  narrowband  ambiguity  function  (n  s  m)  is  the  characteristic  function  of  the 
Wigner  distribution.  The  "generalized"  ambiguity  function  is  the  characteristic  function 
of  Ae  associated  time-frequency  distribution  Reference  18).] 

If  the  echo  signal  is  from  a  collection  of  moving  scatterers  with  reflectivity 
density  piT,D)  defined  in  such  a  way  that 


*.(<)  =  j  j  P(T.D)^ w,(D(,*T))iL^  . 


(5.8) 


then  it  is  easy  to  show  that 

k./.X'.o) = J  J  p(r.D)  A_fo(r  -f).f  • 

0  -«•  ^  u  /  U 


(5.9) 


(The  measure  dTdD ID  has  important  consequences  later  on  and  is  motivated  in  the 
Appendix.) 
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The  quantity  {e,,fjj  is  known  as  the  reflectivity  image  and  can  be  seen  to  be  die 
actual  reflectivity  density  filtered  by  the  cross-ambiguity  function.  In  practice,  this  image 
is  obtained  by  cross-correlating  the  return  si^al  with  a  family  of  time-scaled  and  shift^ 
test  pulses  and  will  be  examin^  in  more  det^  below. 

The  reflectivity  image  is  generally  coii^lex-valued.  For  many  purposes,  we  are 
more  interested  in  the  intensity  image: 


=  j  J  j  ]p(T,D)p’(T'.D’)y 
0->0-> 


dT'dD'dTdD 
D'  D 


(5.10) 


A  common  model  for  the  reflectivity  density  assumes  point  sources  so  that 


p{T,D)p*{r,D')^a^{T,D)S{T-r,D-D') ,  (5.11) 


where  is  the  target  cross-section  density.  Substitution  of  Equation  5.11  into  5.10 
yields  the  (model-based)  result 


^  i-M  » I  j  <j’(r,o) 


ardP 


(5.12) 


Note  that  failure  of  this  model  (i.e.,  if  the  sources  are  correlated  in  some  sense)  generally 
will  lead  to  image  artifacts  in  the  form  of  cross  terms. 

Equation  5.12  may  be  interpreted  as  a  (perhaps  more  meaningful)  relationship 
between  cross-ambiguity  functions  and  densities.  There  is  a  similarity  between  and 
the  short-time  Fourier  transform.  In  Equation  5.6  the  part  of  the  "window  function"  is 
played  by  ^„{t)  =  s„{t)e‘*“*\  and  so  the  density  appears  as  a  generalization  of  the 
spectrogram: 
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[A^(r,a)f  =aj  'jwl(a{u+t))w^{a{u+t+r))s^{u)sl,{u+T)e~'"*^dudt .  (5.13) 
Comparing  this  with 

j  ^  J  J  J^(e, t)w*(v)k<v  +  dvdddz  (5.14) 

yields  (after  a  suitable  change  of  variables) 

j  ♦(«,  ^  (5.15) 

or 

Jir^(«)s*(«+  x)e~‘*‘du  .  (5.16) 


This  is  the  kernel  for  the  spectrogram  (properiy  time-scaled  as  in  Table  4.1)  and  so  it  is 
instructive  to  compare  Equation  5.15  with  the  kernel  (Equation  4.8)  obtained  for  the 
spectrogram. 

As  discussed  in  Example  4.4,  one  of  the  problems  with  the  spec^gram  is  its 
dependence  upon  a  predetermined  size  of  the  window  function  for  estimating  properties 
of  the  signal.  If  the  signal  varies  slowly  in  comparison  to  the  extent  of  the  w^ow,  then 
the  estimate  will  generally  be  good.  However,  if  the  signal  varies  rapidly  within  the 
window,  then  the  estimate  may  be  poor.  Equation  5.15  shows  that  the  so-called 
scaleogram  of  Equation  5.13  offers  an  interesting  "solution"  to  this  problem.  Here,  a 

fixed-shape  window  s^{t)  wUl  be  time-scaled  by  the  factor  in  such  a  way  that  when 
the  signal  varies  rapi^y  (scale  decreases)  the  window  extent  will  be  correspondingly 
decrea^.  Similarly,  when  the  signal  varies  slowly,  the  extent  of  the  window  will  be 
correspondingly  increased.  This  idi^  is  at  the  heart  of  scale-based  signal  processing  and 
has  bMn  the  focus  of  considerable  research  activity  (References  19  dirough  21). 
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6.  CHOICE  OF  SIGNALS 

Because  of  the  imaging  kernels  ( or  |A„f ),  tl^  main  problem  of  high-resolution 
radar  imagery  is  that  of  choosing  a  pulse  wjt)  and  a  radar  signal  wjt)  that  are  most 
iq>propriate  to  a  specific  imaging  application. 

Recall  that  we  are  concerned  with  die  envelope  s(t)  of  the  waveform  w(f)  =  s(t)e‘"*'. 
For  actual  radar  systems,  both  the  pulse  and  the  signal  will  be  real-valued  and  we  can 
think  of  (t)  and  s^(t)  as  being  the  complex-valued  waveforms  whose  real  parts  are  the 
actual  pulse  envelope  h^{t)  and  radar  signal  envelope  A, (r),  respectively.  For  such  h{t), 

we  have  h{~o>)  =  h*(o))  so  that  h(t)  is  completely  determined  by  its  positive  spectrum. 
It  is  useful  to  define 

s{t)  s  2 J  dot ,  (6.1) 


where  the  real  and  imaginary  parts  of  5(r)  are  a  Hilbert  transform  pair  and  the  factor  of  2 
guarantees  that  s(r)  and  h{t)  have  the  same  energy.  [When  the  spectrum  of  s{t)  has 
positive  support,  then  the  quadrature  model  for  w(f)  yields  the  same  results 
(Reference  22).]  Consequently,  we  can  identify  the  Hardy  space 

H^{R)  =  {se  L^(R):supp(s)  c  [O.oe)}  with  the  space  of  all  real-valued  functions  in  l}{R) 
(Reference  19).  H^{R)  is  a  closed  subspace  of  L^(R). 

The  ideal  imaging  kernel  will  be  one  that  is  sharply  peaked  at  (T,a)  =  (0,1)  and  zero 
everywhere  else.  There  are  limits  to  how  well  this  may  be  realiz^,  however,  and  it  is 
well  known  that  matched  Gaussians  are  the  ’’optinud”  w^{t)  and  w^{t)  for  narrowband 
ambiguity  functions  (optimal  in  the  sense  that  diey  minimize  the  uncertainty  relation). 
Unfortunately,  Gaussians  do  not  belong  to  H^{R).  What  is  required  for  any  good  choice 
of  w{t)  is  (Reference  19): 


i)  wgH\R) 

ii)  w{t)  satisfies  the  admissibility  condition  [1^— ^ 

i  to 


■doxee.  [As  we  shall  see 


below,  this  is  required  if  we  are  to  recover  p{T,D)  from  {e,,f„){t,a),  or  <^{T,D)  from 
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Note  that  while  a  Gausnan  does  not  belong  to  the  function  e'**'  aq)  (-ir*)  is 

very  close  to  being  admissible  when  (o„  is  sufficiently  large.  This  suggests  that  matched 
Gaussians  (or  matched  Gaussian-like  functions)  can  be  expected  to  hold  an  analogous 
role  in  the  theory  of  wideband  ambiguiQ'  functions.  In  fact,  Grossman  and  Morlet 

(Reference  19)  have  identified  the  envelope  that  ”...  plays — ^in  the  //^-theory  that  we  are 
concerned  with — ^the  same  privileged  role  that  the  Gaussian  plays  in  -theory."  For 
>  0,  the  function  is  defined  by  its  Fourier  transform  as 


^  0 


forfi)>0, 
for  nr  ^0. 


(6.2) 


[Cp  is  a  normalization  constant:  Cj  =^|pJxeJq)(-‘l/4fi).] 

If  we  choose  s^(t)  and  »v,(/)  =  e'***'j,(r),  then  the  ambiguity 

function  for  these  waveforms  becomes 


A^iz.a)  Va  J  e*™"  exp(-{P^  In^  a>  -  ln^(a(o))do> .  (6.3) 

0 


Expanding  the  ln^{a(o)  term  in  the  argument  shows  that  the  waveforms  will  be 
matched  if  In  this  case,  we  obtain  some  simplification  and  can  write 


(6.4) 


(which  must  be  determined  numerically.)  A  plot  of  |A^(T,a)|^  for  ^  =  1  is  given  in 
Figure  6.1. 
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FIGURE  6.1.  Wideband  Ambiguity  Function  for  Sp{a>)  =  Cpexp[~-\pln^  to) 
Withp-l. 


The  parameter  effectively  determines  the  signal  bandwidth.  This,  of  course, 
determines  fte  resolution  to  which  the  imaging  kernel  can  represent  the  image,  as  can  be 
seen  by  examining  Figures  6.2  and  6.3.  Smaller  P  means  (i^tically)  larger  bandwidth 
with  corresponding  g*%ater  time  resolution  (but  at  the  cost  of  scale  resolution).  For 

comparison  purposes,  the  narrowband  ambiguity  function  for  Sp{e>)  =  is 

displayed  in  Figures  6.4  through  6.6. 
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FIGURE  6.2.  Wideband  Ambiguity  Function  for  5^(o))  =  Cpexp{^^pln^  to) 
Withes  10. 
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FIGURE  6.3.  Wideband  AmbiguiQr  Function  for  aq?{-^pin^  o) 

With  ^  =  1/10. 
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FIGURE  6,4.  Narrowband  Ambiguity  Function  for 
Sp(Q})  =  D^ejqj(-^fio>^)  With  ^  =  1. 
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FIGURE  6.6.  Narrowband  Ambigui^  Function  for 
exp(~jfia)^)  With  - 1/10. 


Of  course,  to  be  accurate,  we  really  want  neidier  the  reflectiviQr  image  {e,,fj){t,a) 

nor  the  intensity  image  rather,  we  desire  the  reflectivity  density  p{T,D)  or  the 

cross-section  density  a^{T.D).  These  densities  are  directly  related  to  the  target 
properties  alone,  and  to  get  at  Aem  we  must  choose  signals  that  lead  to  imaging  kernels 
that  are  invertible. 

The  reflectivity  image  is  given  by 


(6J) 
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Qeariy,  we  need  to  find  some  such  that 


(6.6) 


If  we  write  B^{r,  a)  in  the  form 


Bjtit.a)  =  Va  J Wy(a(M+  x))wl{u)du  , 


(6.7) 


then  it  is  straightforward  to  show  that  our  inversion  requirement  becomes 


p(r',D')  I  ]p(T.D)\(d'{T  - 


\drciD 

)  D 

(6.8) 


If  =w.,  then  the  admissibility  condition  shows  that  this  is  a  well-defined 
operation.  However,  we  are  back  to  the  original  problem  of  seeking  an  ideal  ambiguity 
function  A^{x,a)  that  is  sharply  peaked  at  (T,fl)  =  (0,1)  and  zero  everywhere  else. 


7.  WHY  GO  WIDEBAND? 


While  the  constraint  on  the  duration-bandwidth  product  (the  uncertainty  principle) 
guarantees  that  we  can  never  achieve  the  optimal  imaging  kernel  \{x,a)  =  S(T)S(a  - 1), 

we  are  able  to  obtain  A^(T,a)  =  S(t}  or  A^{x,a)  =  5{a  - 1).  In  fact,  since  we  can  make 
multiple  radar  measurements  of  the  same  system,  we  could  follow  a  measurement  using 
the  waveform  described  in  Figure  6.2  (Figure  6.5)  with  a  second  measurement  using  the 
waveform  of  Figure  6.3  (Figure  6.6).  Correlating  such  ideal  measurements  allows  us  to 
achieve  as  much  imaging  resolution  as  we  might  desire.  In  either  of  these  measurements 

we  are  free  to  choose  from  a  wide  variety  of  wavefomns  w(/),  since  we  no  longer  seek  to 
"optimize"  against  some  uncertainty  constraint  This  freedom  in  choice  of  w{t)  means 
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that  we  can  select  a  wavefoim  that  is  suited  to  other  aspects  of  the  particular  problem  at 
hand. 

Even  when  we  are  not  planning  such  correlated  measurement  imaging,  we  still  are 
particularly  interested  in  high-resolution  downrange  imaging,  since  this  information 
might  be  used  more  readily  for  target  imaging  tiian  cross-range  (Doppler)  information,  hi 

either  case,  however,  we  are  interested  in  those  w(t)  that  have  very  narrow  envelopes 

s(r),  since  these  so-called  wideband  waveforms  allow  for  very  fine  downrange  resolution 
of  radar  targets.  (In  contrast,  narrowband  signals  provide  low  downrange  resolution.) 
Much  of  the  intuition  developed  fiom  experience  with  ordinary  (narrowband)  radar 
systems  carries  over  directly  to  the  wideband  case.  In  general,  however,  the  narrowband 
amlnguity  function  may  not  be  a  completely  accurate  imaging  kernel  for  rapidly  moving 
targets.  Note  also  that  the  narrowband  ambiguity  function  is  related  to  the  Wigno* 
distribution,  while  the  wideband  ambiguity  function  does  not  enjoy  the  same  relationship. 
However,  we  can  dismiss  any  vague  concern  that  Fourier  methods  may  not  apply  to 
wideband  radars.  As  we  have  seen,  Fourier  methods  work  just  fine  in  the  present 
analysis — provided  that  they  are  applied  correctly. 

To  understand  better  the  differences  between  the  two  approaches  and  to  quantify  the 
possible  loss  in  (narrowband)  imaging  accuracy,  we  shall  compare  A^{x,a)  with 

various  bandwidths.  If  we  set  w,(/)  **  H'^(r)  =  e‘***'s{t)  then,  in  terms  of  the 
narrowband  imaging  kernel,  we  can  write 


X(t.{a  -  l)a)J  s  2;rVa  e"''  xjt,{a  -  IK) 


=  e-''  J s{u + r)s*(ii)/‘"')**“  du 

=  Vac'**'  J s(©)5*(© - (1  - a)fl),)e'* do). 


Whereas,  the  corresponding  wideband  imaging  kernel  is 


A(f  la,a)s  A^(t  la,a)-^ e'**'  J s{au + du 


(7.2) 


=  Vo  e'*''  J s(ca)V{aa)  -  (1  -  a)(0, )e**  don. 
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As  a  concrete  example,  consider  the  Gaussian  waveform  s(a})  =  Dpexp[-^P(0^). 
This  waveform  is  computationally  easy  to  deal  with  and  is  almost  in  H^(R)  when  is 
sufficiently  large — say  eo,  ^  Bandwidth /2.  [In  this  case,  it  is  appropriate  to  use  the 
Woodward  Bandwidth  (Reference  23): 


«•  «• 

=  J dm  GJ^e 


2)5 


(7.3) 


since =‘s[fi71c.]  Substitution  of  this  Gaussian  waveform  into  Equations  7.1  and  7.2 
yields 


x{t,(a-l)to„)  =  ^e‘“‘‘exp(~^p(l-afml)]exF{-pa)^+{pG)„i\ -a)+it)o>]dQ} 


-  ^fa  ^  <u,r jctp(-^^(l  -  af  ml  )exp  -- 


t^-{\-afp^m, 


> 


(7.4) 


and 


exp(-^Pi^  -  afo>l)  J  «p[-i)5(l + a^)m^  + 


{pm„a{l  +  a) + it)m^dm 


(7.5) 

r^-a^(l-a)^^V' 

2)3(1  + a*)  / 


For  fixed  a,  the  scatterer  location  (in  tuiK)  will  be  estimated  by  the  maximum  of  the 
respective  ambiguity  functions.  If  denotes  the  estinuite  from  the  wideband  ambiguity 

function  A{t I a,a)  and  the  narrowband  ambiguity  function  x{t,{}.-a)m^)  estimate, 
then  it  is  easy  to  see  that 
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=  (7.6) 

for  0),  =  \B„ .  This  difference  vanishes  when  the  target  is  not  moving  (fl  =  1).  However, 
downrange  distances  of  moving  targets  will  be  overestimated  by  the  narrowband 
ambiguity  function.  But  this  error  is  usually  not  significant  if  the  target  velocity  obeys 
|v|«c;  hence. 


(a 


4B^c^  ‘ 


(7.7) 


(We  can  show  that  the  phase  error  is  similarly  insignificant,  but  the  phase  term  is  usually 
ignored.) 

In  the  narrowband  case,  it  is  stiU  true  that  increasing  bandwidth  will  increase  time 
resolution.  The  example  shows  that  when  s(a))  is  chosen  to  be  symmetric  and  centered 
on  a)„  ^  there  is  little  difference  between  the  wideband  and  narrowband  results  (at 
least  as  far  as  radar  systems  are  concerned).  This  is,  perhaps,  why  the  "widebandness"  of 
radar  is  sometimes  defined  by  the  ratio  of  bandwidth  to  center  frequency. 

In  the  Hardy  space  case,  no  such  minimum  size  restrictions  for  a,  exist;  however,  it 
is  important  to  remember  that  these  waveforms  should  be  elements  of  Hardy  space.  TUs 
requirement  guarantees  that  we  vdll  not  have  any  unexpected  artifacts  obtained  by  trying 

to  represent  functions  that  are  members  of  H^{R)  by  elements  of  L*(/?).  Moreover, 
often  it  may  be  advantageous  to  consider  iH’oper  waveforms  (i.e.,  with  nonsymmetric 
spectra)  for  radar-detection  applications  (c.f..  Reference  21).  In  this  case,  the  differences 
between  the  narrowband  and  wideband  imaging  kernels  may  become  significant  because 
of  the  way  the  imaging  kernel  is  sharpened  when  a  <  1 . 


8.  CONCLUSION 


It  is  not  surprising  that,  for  radar  applications,  the  wideband  and  narrowband 
ambiguity  functions  should  yield  such  close  results.  After  all,  they  only  differ  by  terms 
of  order  vie.  Moreover,  if  there  were  serious  errors  in  the  way  the  namwband 
ambiguity  functions  performed,  then  surely  they  would  have  been  discovered  in  a  half- 
century  of  radar-system  development  (Sonar  applications  are  another  matter,  however. 
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and  all  of  these  results  cany  over  to  acoustic  sensing  but  with  the  possibility  that  vie 
can,  in  practice,  become  significantly  large.) 

However,  we  should  not  conclude  that  wideband  ambiguity  functions  are  merely  a 
mathematical  curiosity.  The  narrowband  analysis  does  not  tell  the  entire  story,  and 
intuition  gained  from  it  is  not  complete.  The  shift  in  attention  from  time-frequency  to 
time-scale  representations  is  made  natural  by  the  shift  from  narrowband  to  wideband 
ambiguity  functions.  In  turn,  the  representation  of  sisals  by  affine  covariant  states 
(wavelets)  holds  the  potential  for  si^fficandy  improving  the  imaging  capabilities  of 
future  rad^  (and  sonar)  systems  (c.f..  Reference  21  and  references  cited  therein). 

Current  work  has  been  concentrated  on  demonstrating  the  utility  of  this  approach 
using  actual  and  simulated  frequency  domain  data.  Stepps  fluency  data  collected 
over  large  bandwidths  have  been  in  use  for  many  years,  and  filtering  techniques  allow  for 
easy  mt^fication  to  the  present  purposes. 

Current  plans  for  future  work  include  a  detailed  investigation  into  appropriate 
choices  for  radar  waveforms,  as  well  as  {nactical  techniques  for  generating  them  in  the 
time-domain. 
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Appmdix 

MOTIVATION  FOR  THE  MEASURE  dTdD  /  D 


In  the  nanowband  case,  the  equivalent  tt)  Equation  S.9  can  be  writtra 

«•  m 

(e./X/.fl))=  J  \p{T,a)x{T-t,Q-(0)drdn  .  (A.1) 


This  is  a  correlation  integral  between  p  and  so  obeys  the  (Fourier) 

correlation  theorem: 

]gix')h{x'+x)dx'  ^  g{y)h\y).  (A.2) 


where  ^  denotes  Fourier  transformation.  We  would  like  Equation  5.9  to  also  obey  a 

correlation  property.  The  measure  choice  dx/x  aUows  us  to  do  this,  since  under  the 
Mellin  tranrform  (Reference  24), 


8{y)K-y) .  (A-3) 


where  <=>  denotes  the  Mellin  transformation: 

M 
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iiy)  *  / « dx  .  (A.4) 

0 

Altonately,  we  can  avoid  any  mention  of  Mellin  transforms  and  force  integrals  of 
the  type  in  Equation  A.3  to  appear  to  be  the  kind  of  correlation  integrals  we  are  used  to. 

This  can  be  accomplished  by  making  the  change  of  variables  x  e“,  which  yields 


(A.5) 


where  G{u)  s  (This  approach  is  used  in  Reference  19.) 

Independent  of  the  motivation  for  this  choice  of  measure,  the  consequences  have 
been  important  for  us.  It  was  this  measure  that  was  responsible  for  the  signal 
adntissil^ty  condition  of  Section  6  and  for  the  logarithm  spearing  in  the  argument  of 
Equation  6.2. 
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